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A diffusion-induced transition in the phase separation 
of binary fluid mixtures subjected to a temperature ramp 
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Abstract. - Demixing of binary fluids subjected to slow temperature ramps shows repeated waves 
of nucleation which arise as a consequence of the competition between generation of supersatura- 
tion by the temperature ramp and relaxation of supersaturation by diffusive transport and flow. 
Here, we use an advection-reaction-diffusion model to study the oscillations in the weak- and 
strong-diffusion regime. There is a sharp transition between the two regimes, which can only 
be understood based on the probability distribution function of the composition rather than in 
terms of the average composition. We argue that this transition might be responsible for some yet 
unclear features of experiments, like the appearance of secondary oscillations and bimodal droplet 
size distributions. 



Introduction. — Liquid-liquid phase separation oc- 
curs whenever an isotropic binary mixture is transferred 
into a bi-phasic region where its isotropic state is no longer 
stable, and the mixture decomposes into two equilibrium 
phases. In classical approaches [I] demixing was induced 
by sudden temperature quenches, after which the phase 
separation was monitored in constant temperature condi- 
tions. In many industrial [5] and natural [3]|1] applications 
however, temperature is not constant but it coevolves with 
' the phase separation. Accordingly, some recent works fo- 
cused on systems subjected to time-dependent variations 
of the temperature, like repeated cycles of cooling and 
heating [5] , or slow continuous temperature ramps [6Ul3j . 

Experiments [BHH] in fluid mixtures subjected to slow 
temperature ramps show unexpected phenomena: rather 
than continuously following the gradual change in tem- 
perature, the phase separation exhibits consecutive bursts 
of droplet nucleation alternating with quiescent periods 

|5HI3|. 

This intriguing phenomenon can be explained in terms 
of the competition between the temperature ramp and dif- 
fusion. Due to the temperature ramp there are droplets 
present in the system at any time. Any change of temper- 
ature results in a change of the equilibrium composition 
of the background fluid and of the droplets: supersatu- 
ration - defined as the deviation from the actual to the 
equilibrium composition — builds up in the sample. This 



supersaturation is relaxed by diffusive exchange of mass 
between the droplets and their environment, or occasion- 
ally also by nucleation of new droplets. 

Immediately after a wave of nucleation, the droplet con- 
centration is high [8 . In this case, diffusion, which acts on 
length scales of typical droplet distances, relaxes supersat- 
uration by orders of magnitude faster than it is generated 
by the temperature ramp |12j . The concentration remains 
close to equilibrium, and no new droplets are nucleated. 
The system enters a quiescent period. In the course of 
time, the number of droplets in the system decays due 
to coarsening of the droplet distribution and sedimenta- 
tion. Droplet distances grow. Thus the diffusive exchange 
becomes slow, and supersaturation rises again. Eventu- 
ally, a new wave of nucleation is triggered. Clearly the 
strength of diffusive transport has a critical impact on the 
alternation of the active and quiescent periods of phase 
separation. 

The strength of the diffusive transport is characterized 
by the ratio D/L 2 . Hence, it is affected not only by the 
distance between droplets L, but also by the diffusion coef- 
ficient D. The latter effect is especially important because 
the temperature ramp, that drives the phase separation, 
causes considerable variations of D. It rises monotonously 
as function of the distance from the critical point (cf. [14]). 

Here, our aim is to investigate systematically the be- 
haviour of the oscillations in different (weak and strong) 
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diffusion regimes. We adopt a numerical approach built 
on previous mean-field discussions [HUH] as iar as the evo- 
lution of the supersaturation is concerned, and combined 
with a reactive-flow description [T3]. The latter deals ex- 
plicitly with the spatial distribution of droplets and super- 
saturation, and with the dynamics of droplets (formation, 
coagulation and advection). Upon increasing the diffu- 
sion coefficient we find a clearly defined transition in the 
characteristics of the oscillations — in particular in their 
period. We will compare this behaviour to experimental 
results, and point out that some yet unclear features of the 
experiments might be attributed to the continuous growth 
of the diffusion coefficient as the system departs from the 
critical point. 

The model. — It is convenient to describe the dynam- 
ics in terms of a normalized composition a |11H13| , defined 
in such a way that it takes the value Co = 1 in equilib- 
rium, and values < a < 1 when the system departs from 
equilibrium. In experiments, the normalized composition 
decreases like a = — £<r, where the decay rate £ is a func- 
tion of temperature [7]fT0] . It can be kept almost constant, 
however, by choosing a temperature protocol that assures 
that the equilibrium flux onto the droplets is constant in 
time Q3QI]. 

In order to focus on the spatial fluctuations of the nor- 
malized composition we formulate the demixing dynamics 
in terms of the scalar field o~(x, t). It is subjected to advec- 
tion by the background flow, to an exponentially decaying 
reaction due to the temperature ramp, and to diffusion 
(for details see |13j): 



da(x,t) 
dt 



-u-Vcr(x,t) = T(o)-£a{x,t)+VAa{x,t). (1) 



Here V is the diffusion coefficient, and u the velocity field 
of the background convection. The term T(o~) represents 
formally the nucleation of droplets (see Eq. (TJJ for the 
actual implementation as a stochastic process). Nucle- 
ation brings the system back to equilibrium by resetting 
the local composition to do = 1. Thus it brings about a 
non-trivial (i.e. non-zero) final distribution of the o~{x, t) 
field. 

For numerical investigations of the advection-reaction- 
diffusion equation, Eq. (JTJ, we use a lattice model pro- 
posed and described in an earlier work [13j . Here we just 
outline the basic elements of the model. 

Spatial degrees of freedom are coarse grained in the form 
of cells of size e which are labelled by the index (i,j) on an 
N x N lattice defined in a square of size L with periodic 
boundary conditions. At each time r, we apply a reaction 
step in each lattice cell (i,j), 



T (iJ) -> o-{i,j) -£,t a(i,j), 



(2) 



that describes the changes of the normalized composition 
due to the external cooling. This step is followed by a 
diffusion step according to the transformation (see |15j): 



where D(i,j) = -a(i,j) + [a(i + l ) j) + o-(i-l,j) + a(i,j + 
1) +a(i,j — l)]/4 is the discrete Laplace operator, and D 
the diffusion coefficient in the discrete map. In the time 
continuous limit (r — > 0), this D corresponds to a diffusion 
coefficient V = De 2 /A. 

In the next step, we nucleate droplets in cells that are far 
from equilibrium, (i.e. a < o~ t h = 2/3) with a piecewise- 
linear probability a(o~) of nucleation (as in |13j): 



a(a) = 




if a > ath = 2/3 
if 1/3 < a < 2/3 



(4) 



if a < a* 



1/3. 



This choice reflects on one hand that for o < o~ sp spinodal 
decomposition causes instantaneous phase separation [I] 
IT2"] , and on other hand that above a threshold value ath, 
the mixture is so close to equilibrium that practically no 
nucleation events occur. Upon nucleation the composition 
is set back to ao(x,y) = 1 in each cell where nucleation 
took place, as well as in the eight neighbouring cells. 

The effects of the background convection are repre- 
sented by an advection step that mixes the lattice cells 
and droplets in the lattice. This is generated by an incom- 
pressible alternating shear flow (as in [16]) that produces 
chaotic trajectories. The time-periodic velocity field of the 
flow is given by 



u x {x,y,t) = 4± 9(| -imodT) sin 2 ^, 
u y (x,y,t) 



4± e(tmodT-f) sin 2 ^ 



(5) 



o-(i,j) -+a(i,j) +Dt- D(i,j), 



(3) 



where Q(x) is the Heaviside step function, T the period of 
the flow, and A the strength of the flow. 

Finally, we take into account the merging of droplets: 
whenever two droplets approach each other to within a 
distance r$ , the two droplets are replaced by a new droplet 
placed in the centre of mass of the two droplets. 

Henceforth, we fix the cell size to the specific value e* 
that represents the maximal spatial scale where the a- 
field takes constant values in the presence of diffusion and 
flow. Below the scale e* the small-scale fluctuations in a, 
that are created continuously by the straining action of 
the velocity field, are dissipated by diffusion [T7]. For a 
flow with Lyapunov exponent A the scale e* can be esti- 
mated as follows |T7|. The stretching and folding caused 
by the chaotic flow generates inhomogeneities of the a- 
field on scales that decrease according to e = — Xe. On the 
other hand, diffusion causes a broadening of the homoge- 
neous areas according to i = D/e (adopting the solution 
e ~ t 1 / 2 ). Thus the size of homogeneous areas evolves 
as: e = D/e — Xe. This equation possesses a steady state 
solution e* = {D/X) 1 / 2 = LPe^ 1 / 2 , where Pe = XL 2 / D 
is the Peclet number, a measure of the strength of the 
diffusivity. Thus, the dimensionless cell size amounts to 
e* = LPc~ ' . In what follows, we set the length and time 
scales to L = 1 and T = 1. 

Oscillations in the weak- and strong-diffusion 
regimes. — To explore how the increasing diffusion co- 
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Fig. 1: Evolution of the droplet density rid in numerical simulations for flow amplitude A = 0.8 and cooling rate £ = 0.04 
in the (a) weak (D — 0.05 - red, D — 0.1 - green line), (b) intermediate (D = 0.14,0.15 and 0.16 in decreasing order of the 
amplitudes), and (c) strong (D — 0.25, 0.35, 0.45, and 0.55 in decreasing order of the periods) diffusion regimes. All simulations 
were started with a uniform initial composition field <to = 0.667. The black thin line represents the diffusionless case. For 
intermediate values of D, as shown in panel (b), secondary oscillations appear. 



efficient modifies the oscillations of the droplet density 
we fix the flow rate, the decay rate of the composition 
and the droplet interaction radius to A = 0.8, £ = 0.04, 
and ro = £*, respectively, and study the oscillations for 
different diffusion coefficients ranging from D — 0.05 to 
D = 0.55. We follow 100 individual realizations of the 
system and average the results. 

The oscillations are qualitatively different for low, inter- 
mediate and high values of the diffusion coefficient. In the 
weak-diffusion regime, Fig. QJa), tne period of the oscil- 
lations increases with increasing D, while their amplitude 
decreases faster than in the diffusionless case. As the diffu- 
sion coefficient increases further, this tendency stops at a 
certain value of D. For intermediate values of D [Fig.QJb)] 
secondary oscillations appear. Subsequently, the period of 
the oscillations becomes smaller again [Fig.QJc)], and the 
decay rate of the amplitude decreases with increasing D. 

Transition between the weak- and strong- 
diffusion regimes. — The period of the oscillations 
takes its maximum value approximately in the middle 
of the range where secondary oscillations are observed 
[Fig. [2a)]. Furthermore, the appearance of secondary os- 
cillations in the intermediate diffusion regime is accompa- 
nied also by dramatic quantitative changes of the decay 
rate of the oscillation amplitude and of the decay of the 
composition in the quiescent period of the oscillations. 

We quantify the decay of the oscillation amplitude by 
following how the difference between the maximum droplet 



density and its average value Arid = n<i,j 



■ rid av evolves 



in time, our data show that for each diffusion coefficient 
this amplitude exhibits an exponential decay with rate 

7amplit i 

Arid ~ exp [-7am P iit(-D) • t] ■ (6) 

As a function of the diffusion coefficient D [Fig.^b)], this 
decay rate has an extremum localized at the end of the dif- 
fusion range in which secondary oscillations are present. 
The fastest decay of the amplitude corresponds to that 
value of the diffusion coefficient at which secondary oscil- 



lations disappear from the system. For higher diffusion 
coefficients, the exponent 7 am piit decreases slightly. Then 
it seems to saturate. 

The most abrupt change can be observed in the be- 
haviour of the average composition in the system. We 
observe that in the quiescent period of each oscillation the 
average composition still (cf. [13]) decreases exponentially 
as 

cr av ~exp[-7 CT (L>) -t] , (7) 

due to the change of temperature. In the absence of 
droplets this decay is determined by the local decay, 
a(x) ~ exp(— £ • i), of the composition in individual cells, 
such that 7 CT = £. When droplets are present in the sys- 
tem the diffusion of supersaturation into the droplets slows 
down the decay of the average composition. For small val- 
ues of the diffusion coefficient the magnitude of this effect 
is small, and the average composition still decays with a 
rate j a ~ £. However, when the diffusivity is increased to 
values where secondary oscillation appear, the decay rate 
of the average composition abruptly drops to a much lower 
value [Fig.E^c)]. 

We conclude that the region where secondary oscilla- 
tions are present represents a transition regime between 
the weak- and strong-diffusion regimes, where the tran- 
sition point D tr is defined by the maximal value of the 
oscillation period. To determine the dependence of the 
transition point on the different parameters of the system, 
we run numerical simulations for different decay rates £ of 
the composition, and for different flow rates A. The po- 
sition of the transition point is not affected significantly 
by the flow, but it is very sensitive to the decay rate £. 
More precisely, the transition occurs around D tr ~ 0.16 
for decay rate £ = 0.04, at D tr « 0.08 for f = 0.02, and 
around D tr ~ 0.04 for £ = 0.01. This is compatible with 
a linear scaling 

Ar - 4£. (8) 

The diffusion coefficient D and the decay rate £ define the 
lenght scale A = y/Dj\ at which the effects of the tem- 
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Fig. 2: Diffusion induced transition in the oscillatory behaviour of the phase separating mixture. For £ = 0.04 and A — 0.8 the 
transition occurs in the highlighted region: 0.1 < D < 0.2, determined as the region where secondary oscillations are observable. 
In the transition region: (a) the oscillations' period reaches its maximal value, (b) the decay rate 7 am piit of the amplitude reaches 
its maximum, and (c) the decay rate y a of the average composition er a v decreases suddenly from j a ~ £ = 0.04 to "/& pb 0.03. 



pcraturc ramp and the effects of diffusion are comparable. 
Diffusion can efficiently relax the supersaturation accumu- 
lated by the continuous temperature ramp up to a length 
scale of order A. In the continuum limit we find for the 
transition from the weak to the strong diffusion 



Atr = 




e = e 



(9) 



Hence, the transition occurs when the range in which dif- 
fusion can relax the accumulated supersaturation equals 
the length scale e* = L/Pe ' where the effects of mix- 
ing by advection and diffusion are of the same order of 
magnitude. 

Distributions of supersaturation. At the onset 
of the nucleation wave the average composition a av varies 
in a broad range of values, and does not always approach 
the nucleation threshold a ^ a t h [see Fig.[3]panels (1), (2) 
and (3)]. This indicates that the onset of the nucleation 
wave is not determined by the temporal evolution of the 
average composition. Rather, the nucleation wave is ini- 
tiated by a small fraction of high-supersaturation "spots" 
that reach the threshold. Hence, an adequate description 
of the phenomenon should deal with the full distribution of 
the composition rather than only with the average com- 
position. In this section, we present details about the 
behaviour of the composition distributions. 

We focus on the evolution of the probability distribu- 
tion functions (pdf-s) of the composition field in the low, 
intermediate- and strong-diffusion regimes. To obtain the 
pdf-s we consider composition bins of size 6a = 1/200. 
The resulting pdf-s are presented in Fig. [3] for different 
time instances: the upper panels (1) represent the weak- 
diffusion regime D = 0.05, while the middle (2) and bot- 
tom panels (3), the intermediate D = 0.15, and strong- 
diffusion D = 0.55 regimes, respectively. 

The initial condition is a uniform composition field 
a(x, 0) = 0.667, that is represented in the pdf-s as a 6- 
function (not shown). Due to the temperature ramp, this 



peak moves to the left, in the direction of the nucleation 
threshold ath = 2/3, while it takes the form of a Gaussian 
distribution. When few points of the distribution cross 
the threshold level the nucleation wave starts. The points 
(cells) in which nucleation took place move up to a = 1. 
In the panels (la), (2a), and (3a) of Fig. [3] we show the 
distributions at the end of the first nucleation wave. Re- 
gardless of the diffusion regime the behaviour of the cr-field 
is similar in all the three cases: after several cells nucle- 
ate droplets the distribution becomes bimodal. There is 
a big number of cells being already around the equilib- 
rium value <7 = 1, and a small number of cells remain still 
slightly above the threshold and form a second peak. 

This small peak is subjected to two opposing forces. Dif- 
fusion, represented by the red arrows (ii) in Fig. [3l tends 
to move the peak to the right: i.e., into the direction of 
the average composition that is already close to equilib- 
rium. The second force is that of the temperature ramp. 
It is represented by the black arrows (i) in Fig. |3l and 
favours nucleation, i.e., it drives the peak further to the 
left in the direction of the nucleation threshold. The com- 
petition between the effects of diffusion and temperature 
ramp determines the fate of the small peak. Accordingly, 
the character of the oscillations differs in the three diffu- 
sion regimes. 

In the weak-diffusion regime the effects of diffusion are 
small [(ii)<(i) in Fig. [3] panel (la)]. Thus the tempera- 
ture ramp will dominate, and the peak will slowly cross 
the threshold line a t h and nucleate droplets. This delayed 
nucleation results in a long tail of the main peak, modi- 
fying its Gaussian shape [see Fig [3] panel (lb)]. The fol- 
lowing oscillations [Fig. [3] panels (lc) and (Id)] will take 
place according to the same scenario. The slowing down 
is more and more efficient as the diffusion coefficient in- 
creases. Thus, the duration of the nucleation wave in- 
creases, explaining the increase of the oscillation period in 
the weak diffusion regime. 

In the intermediate diffusion range the effects of the 
temperature drift and of diffusion balance each other 
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Fig. 3: Evolution of the average composition <r av (left panels), and probability distribution functions of the composition at time 
instances marked by points in the left panels in the weak- (D — 0.05 - upper panels), intermediate- (D = 0.15 - middle panels), 
and strong-diffusion regimes (D — 0.55 - bottom panels). The flow rate is A = 0.8, the decay rate £ = 0.04. The simulations 
were started with a uniform initial composition field <tq = 0.667. The plots of the pdf-s all show the range 0.5 < cr < 1, the 
arrows (i) and (ii) are described in the text. 



[(ii)~(i) in Fig. |3] panel (2a)]. The small peak is arrested 
from nucleation and stopped in the close vicinity of the 
nucleation threshold without actually reaching it. Mean- 
while, the continuous ramp of temperature moves the big 
peak to lower composition levels where the effect of dif- 
fusion becomes less efficient. Eventually, the small peak 
manages to cross the nucleation level. This happens after 
the first nucleation wave was finished and gives rise to a 
separate small wave of nucleation. In this way, a secondary 
oscillation of the droplet density [as shown in Fig.QJb)] oc- 
curs in the system. After this secondary nucleation wave a 
distinct small peak forms close to a — 1 which does not at- 
tach to the big peak. Both peaks maintain their Gaussian 
character while travelling downward [Fig. [3] panel (2b)], 
maintaining the bimodal distribution of the composition. 
In the next wave of nucleation, the bigger peak breaks 
again in two distinct modes, forming a trimodal distri- 
bution. The number of modes increases by one at each 
subsequent principal or secondary nucleation wave. Thus, 
the distribution converges rapidly to a flat distributions in 
which the standard deviations of the individual modes are 
larger than the distances between them, and the phase 
separation proceeds continuously rather than in distin- 
guishable oscillations. This explains the rapid decay of 
the amplitude in this transition regime. 

In the strong-diffusion regime the effects of diffusion 
overcome the effects of the temperature ramp [(ii) > (i) in 
Fig.[3]panel (3a)]. The small peak is not only arrested from 



nucleation, but it moves backwards towards the big peak. 
The big peak in turn, travels downward in the direction 
of the nucleation threshold because of the intense diffusive 
exchange. Soon the two peaks merge [Fig. [3] panel (3b)] 
and form again one single peak with an (almost) Gaussian 
shape [Fig. [3]panel (3c)]. In the next nucleation waves ev- 
erything is repeated in a similar way [Fig. [3] panels (3c) 
and (3d)]. In this regime, there is always a certain number 
of cells in which the supersaturation is released by diffu- 
sive transport. For high values of the diffusion coefficient, 
the downward travel of the big peak is accelerated by the 
rapid diffusive exchange with the small droplets. It hence 
moves much faster than in the weak diffusion regime where 
this motion was caused by the temperature ramp alone. 
As a result the period of the oscillations decreases in this 
regime. 

When starting the simulations with a random initial 
condition instead of uniform a er-field this picture does not 
change because diffusion tends to synchronize the compo- 
sition. In the strong diffusion regime the synchronization 
of lattice cells is so fast that the character (period, am- 
plitude) of the oscillations is hardly affected by the initial 
condition. Even in the weak diffusion regime the oscil- 
lations are still clearly observable. Their period remains 
almost the same, while their amplitude is considerably re- 
duced as compared to the case of the uniform initial con- 
dition. 
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Discussion and conclusions. — Oscillatory phase 
separation of binary mixtures subjected to slow tempera- 
ture ramps is a relatively new topic in liquid-liquid phase 
separation. Even though precision experiments unveiled 
many characteristics of the phenomena [BHE], and there 
are numerous theoretical attempts to explain their dy- 
namics [5T-TT3"]. the parameter dependence of the oscilla- 
tions period is still an open question. Secondary effects 
appear even in well-controlled experiments 8, Fig. 15], 
and shield the principal behaviour of the system. For in- 
stance, (a) secondary oscillations can appear [15] . (b) the 
frequency of the oscillations increases during one experi- 
mental run (Fig. 10 in f ; Fig. 3 in QT]; Fig. 15 in [SJ), and 
(c) bimodal droplet size distributions occur [3] Fig. 3]. 

In this communication we identified these secondary ef- 
fects in a model system, and we pointed out that they 
can be attributed to the growing diffusion coefficient as 
the system departs from the critical point. Upon increas- 
ing the diffusion coefficient we observe a cross-over regime 
where the behaviour of the demixing changes qualitatively. 

(a) In the cross-over regime the system performs a sec- 
ondary oscillation between any two major bursts of nucle- 
ation [Fig.QTb)]. Similar secondary oscillations are clearly 
visible in experiments when the system is still close to the 
critical point, and accordingly the diffusion coefficient is 
relatively small [18]. The secondary oscillations are ob- 
servable over several periods. In the experiments they 
disappear eventually, as the diffusion coefficient increases 
due to the change of temperature. 

(b) Our result, Fig. [2ja), shows that in experimental 
runs where the diffusion coefficient increases gradually as 
the system departs from the critical point, one should ex- 
pect that the period of the oscillations might initially in- 
crease slightly and eventually decrease considerably. This 
has indeed been observed in experiments (cf. Fig. 10 in 
0; Fig. 3 in [II]; Fig. 15 in [SJ). 

(c) Exploring the distribution of the supersaturation 
has revealed that the oscillatory phase separation can not 
be described adequately in terms of the average composi- 
tion. The spatial distributions of the composition plays a 
crucial role, and should explicitly be taken into account. 
This is particularly important in the cross-over regime 
where the secondary oscillations are observed. They are 
consequences of the bimodal distributions that temporar- 
ily appear in the evolution of the composition field. Ac- 
cordingly, droplets are formed at two different time in- 
stances during an oscillation. Since droplets start to grow 
immediately after they are nucleated, this difference in 
the nucleation time will lead to bimodal droplet size dis- 
tributions as they have been observed in [21 Fig. 3]. We 
expect that such bimodal droplet size distributions appear 
only in the intermediate and strong-diffusion regimes. An 
experimental test of these predictions is under way |18j . 

In summary, we conclude that identifying secondary ef- 
fects which arise from the interplay of diffusion and the 
non-trivial spatial distribution of supersaturation provides 
an important step to a comprehensive understanding of 



the dynamics of phase separation in the presence of a tem- 
perature ramp. 
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